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1 .  INTRODUCTION 

The  performance  r  naval  radars  is  routinely  checked  by 
measurement  of  vertical  ''overage  diagrams.  This  may  be  accomplished  by 
measuring  the  signal  returned  from  a  target  of  known  radar 
cross-section  flying  >-  constant  aspect  profile,  or  from  a  sphere 
dropped  from  a  known  reight  and  range.  Measurements  carried  under 
'standard'  conditions  'an  be  compared  with  nominal  performance  figures 
as  a  test  of  system  d-  radation. 

It  is  not  'i'fi*'\al,  however,  for  radar  returns  several  decibels 
above  or  below  the  ejected  levels  (or  longer/ shorter  detection  ranges) 
to  be  measured  wi  ,hftut  anomalous  propagation  or  system  degradation. 
Much  of  the  disc*  psney  may  be  explained  in  terms  of  variations  in 
humidity,  sea  st»  e  and  the  use  of  different  target  RCS  fluctuation 
models. 


The  following  sections  outline  the  method  and  theory  behind  two 
computer  programs.  The  first  is  concerned  with  the  evaluation  of 
signal,  clutter  and  noise  levels  under  specified  environmental 
conditions  and  target  profile.  A  second  program  uses  these  data  to 
calculate  pxobabil i ties  of  detection  for  appropriate  scaling 
parameters,  opera  tor/ auto-de tection  and  generalized  chi-squared  toiget 
fluctuation  models,  including  Swerling  and  Weinstock  cases.  The 

interactive  versiov.  of  the  programs  is  written  in  BASIC  langua>  >  and 
requires  32K  bytes  of  RAM,  Execution  times  are  approximately  5  seconds 
per  range  point  on  a  Tektronix  4051  Graphics  System. 

2 .  GENERAL  BACKGROUND 

2.1  The  Radar  Equation 

The  radar  equation  relates  the  returned  power  at  the  antenna  (Pr)  to 
the  transmitted  power  (P^)  and  power  gain  (G)  at  wavelength  (X)  for  a 
target  of  cross-section  (a): 

p  _  PtG*^*F*  1 

r  “  (4n)}  R4  L 
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This  equation  differs  from  the  free  space  equation  by  the 
inclusion  of  an  atmospheric  loss  factor  (L)  and  the  pattern  propagation 
factor  (F)  which  accounts  for  multipath  and  diffraction  effects.  An 
exponential  atmospheric  model  is  assumed,  and  multipath  includes  only 
the  direct  ray  and  one  indirect  ray  which  is  reflected  from  the  sea 
surface.  The  geometry  is  shown  in  figure  1. 

In  a  dynamic  atmosphere  F  and  L  are  non-uniform  in  space  and 
time,  leading  to  additional  fading  effects,  however  the  current  model 
calculates  only  average  values.  The  only  sources  of  signal  fluctuation 
considered  are  the  radar  cross-sections  of  the  target  and  sea  surface. 


For  a  given  probability  of  detection  a,  usually  50  or  95%,  the 
maximum  range  Rq  for  noise-limited  detection  is 


'  PtG*oA*F<  1  0.1S 
.  (4n) jDq  L  P„  J 


where  Pfl  is  the  noise  power  and  D0  is  the  signal-to-noise  ratio  required 
to  yield  the  desired  detection  and  false  alarm  probabilities.  The  program 
assumes  Rayleigh  amplitude  statistics  for  clutter  and  receiver  noise,  so 
that  Pn  above  can  be  read  as  noise-plus-cluttcr  power.  Mean  clutter 
power  is  based  on  a  fit  of  Nathanson's  data  (ref  6),  and  the  mean  noise 
level  is  as  described  below. 


2.2  System  Noise  Power 

Equations  1  and  2  are  expressed  in  terms  of  powers  referred  to 
the  antenna,  the  noise  power  Pn  is  also  described  by  a  system 

temperature  Tg  referred  to  the  antenna 

Pn  =  kTsBn  -  kV*  3. 

where  k  is  Boltzmann's  constant,  Bn  is  the  noise  bandwidth,  t 
is  the  effective  pulse  width,  and  the  system  temperature  is 
conveniently  expressed  in  the  form  used  by  Blake  (Reference  1) 

Ts  =  Ta  +  Tr  (Lr  -1)  +  Lr  T0  (NF  -  1) 


4. 
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in  which  T  ,  Lr,  are  the  receiving  line  temperature  and  losses, 
NF  is  the  receiver  noise  figure,  TQ  =  290  and  Ta  is  the  antenna 
temperature.  Ta  is  calculated  by  an  equivalent  path  (Req)  method 

Ta  =  290  +  {  lOVf  +  5.8xl0‘/fJ-s  -  285  )/(Req  y0  )  5. 

where  f  is  in  MHz,  y0  is  the  sea-level  atmospheric  attenuation 
rate  (loss  per  n.mile)  and  R  in  n.mile  is  approximately  related  to 
the  beam  tilt  (p) 

log10(Reo  -  3.5)  =  2.03  +  241ogi0 (1  -  20/n)  6. 

Antenna  temperatures  for  standard  conditions  are  plotted  in 
figure  2  using  equations  5  and  6.  These  approximations,  which  are  not 
described  elsewhere,  agree  to  within  a  few  degrees  with  the  results  of 
a  numerical  integration  described  in  reference  2.  The  atmospheric 
attenuation  is  based  on  the  Van  Vleck-Weisskopf  model  for  oxygen  and 
water. 

2.3  Atmospheric  Losses 

Atmospheric  losses  include  attenuation  by  oxygen,  water  vapour 
and  precipitation.  Molecular  absorption  by  water  and  oxygen  is 
calculated  using  known  intensities  and  line-widths  of  the  rotational 
absorption  lines.  The  relative  contributions  of  water  and  oxygen  are 
shown  in  figure  3.  At  UHF  frequencies  the  absorption  is  dominated  by 
the  tail  of  the  50  GHz  absorption  band  of  oxygen,  which  is  calculated 
using  the  Van  Vleck-Weisskopf  model  (reference  3)  with  semi-empirical 
line-widths  (reference  1)  to  account  for  band  overlap  effects. 

Many  empirical  expressions  and  experimental  data  exist  for 
rainfall  attenuation  above  3  GHz  (references  1-5),  however  this  is 
deliberately  omitted  from  the  program.  At  search  radar  frequencies  the 
attenuation  is  negligible  (10-J  -  10-5  dB/km)  and  the  incorporation  of  this 
loss  may  introduce  errors  larger  than  this  figure  due  to 

(i)  uncertainty  in  the  cut-off  of  the  cloud  base 

(ii)  inhomogeneity  of  the  rainfall  rate  over  the  propagation 
path  and 
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(iii)  complex  frequency-dependence  due  to 

(a)  the  transition  from  Mie  to  Rayleigh  turbidity  with 
variable  drop-size  distribution 

(b)  rapid  dispersion  of  the  dielectric  constant  of  water. 

Volume  backscatter  from  rain  is  more  important  than  attenuation 
at  the  frequencies  of  interest,  but  this  is  simply  another  aspect  of 
the  problem  described  above.  Both  are  treated  in  the  more  general 
radar  model  of  reference  2. 


2.4  Miscellaneous  Losses 

The  probability  of  detection  (Pd)  can  be  calculated  from 
the  ratio  of  signal  energy  to  noise  power  per  unit  bandwidth.  In  the 
program  Pd  is  derived  from  mean  signal-to-noise  power  ratios, 
with  the  assumption  of  approximately  matched  pulse  spectrum/ fil ter 
transfer  function. 

If  the  mismatch  is  small  it  can  be  included  as  a  loss  factor 
which  is  calculated  from  the  empirical  expression 

<t,  ■  “  [&]  •[»  *  fe]1  7- 


where  a  is  the  optimum  time-bandwidth  product  (Bt)  for  the  combination 
of  pulse  and  filter  shape,  and  K  is  the  minimum  loss  which  occurs  at  Bt 
=a.  Using  data  from  Blake  (reference  1)  and  Barton  (reference  4)  with 
equation  7,  bandwidth  mismatch  factors  for  various  pulse/filter 
combinations  have  been  plotted  in  figure  4.  The  default  pulse/filter 
shapes  in  the  program  are  rectangular/rectangular,  and  is  calculated 
for  either  operator-  or  auto-detection. 

A  scanning  lots  factor  il.2  dB)  is  applied  to  the  signal  power, 
but  not  the  clutter.  This  is  because  a  beam-shape  loss  factor  has 
already  been  included  in  the  calculation  of  the  effective  ’imuth 
beamwidth  (0jj  -  0.7S  63^3)  used  to  calculate  the  clutter  cell  area. 


0*1  0-2  0*4  1  2  4  10 

TIME -BANDWIDTH  PRODUCT - *- 

Fig.  4.  Bandwidth  mismatch  loss  (dB)  for  combinations 
of  gaussian  (G)  and  rectangular  (R)  pulse/ 
filter  shapes. 
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The  compressed  pulse  width  of  e  pulse  compession  radar  can  be 
substituted  directly  into  the  radar  range  equations  since  the  time- 
bandwidth  product  equals  the  pulse  compression  ratio  for  a  matched 
filter.  A  loss  must  however  be  included  for  practical  radars  where 
amplitude  tapering  has  been  used  to  reduce  the  first  time/range 
sidelobe  below  the  13  dB  for  a  matched  filter.  The  model  assumes  a  1.5 

dB  processing  loss  based  on  typical  weightings  (e.g.  Hamming  or 

Taylor),  time-bandwidth  products  (Bt  ~  50)  and  mainlobc  pulse  broadening 
(50-60%).  This  loss  is  applied  tc  the  signal  but  not  the  clutter  since 
the  target  is  treated  as  a  point  scatterer  whereas  the  clutter  'fills' 

the  resolution  cell. 

3  PATTERN  PROPAGATION  FACTOR 

3.1  General  Procedure 

The  pattern  propagation  factor  (F)  of  equation  1  is  calculated 
using  (i)  ray  theory  for  the  interference  region,  (ii)  tho  single-mode 

diffraction  result  well  beyond  the  horizon  and  (iii)  an  interpolation 
procedure  near  the  horizon,  A  full  description  of  the  multipath 
calculations  is  given  in  references  1  and  2  with  geometry  and  symbols 
as  per  figure  1.  Annex  III  gives  the  main  geometric  approximations 
used  in  the  BASIC  program.  The  interference  'obing  structure  is  given 
by  equation  1  with  maxima  and  minima  determined  by 

F  =  f (01)  (1  ±  x)  8. 

where 

x  =  9 

fCGj) 

in  which  f(0^)  and  ffGj)  are  the  normalized  antenna  pattern  functions  for 
the  direct  and  indirect  ray  respectively,  r  is  the  surface  roughness 
factor,  ps  is  the  electromagnetic  reflection  coefficient  and  D  is  the 
divergence  factor. 

Beyond  the  horizon  the  single-mode  diffraction  result  is  used  as 
described  in  references  1  and  3.  At  intermediate  ranges  an  exponential 
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interpolation  is  used  to  calculate  F.  The  limits  of  the  intermediate 
region  are: 

Gfflax  =  1.05  x  (radar  horizon) 

Gmin  =  •'  ~  2YC  “  sin-1(ae  cosy^n^  +  ae))  ^ 

sin_1(ae  cosyc/(h2  +  ae))) 

where 

Yc  =  tan“1(X/2nac)J/ J  11. 

Figure  5  is  a  typical  plot  of  F  versus  range  for  a  400  MHz  UIIF 
radar  showing  the  interpolation  region.  The  dotted  line  F=R/Ru  is  the 
single-sc?.n  detection  'threshold'  for  a  target  which  is  known  to  have 
probability  a  of  painting  at  range  Ra  in  the  absence  of  multipath 
effects.  The  single-scan  detection  range  is  found  at  the  intersection 
of  this  line  with  the  curve  for  F. 

A  similar  calculation  for  I-band  is  shown  in  figure  6.  Here  the 
slight  discontinuity  in  the  slope  using  a  simple  interpolation  (linear 
in  dB  versus  range)  is  justified  since  at  higher  frequencies 
clutter-free  detection  is  essentially  horizon-limited.  Over  the 
UIIF/microwavo  region  the  interpolation  procedure  used  should  introduce 
negligible  additional  error  in  the  detection  ranges  compared  with  other 
sources  of  error. 

The  clutter  return  and  signal  level  (1  m1  target)  for  the 
UIIF  radar  and  geometry  described  above  is  plotted  in  figure  7  to 
illustrate  the  significance  of  the  intermediate  region.  At  the  radar 
horizon  the  signal-to-noise  ratio  is  -36  dB  and  increases  by  only  2 
dB/n.mile  until  unity  S/N.  In  this  region  the  probability  of  detection 
is  strongly  dependent  on  target  RCS  fluctuation  -  tho  worst  case  being 
a  steady  target. 

3.2  Reflectivity 

The  depth  of  the  nulls  and  the  enhancement  due  to  multipath  is 
determined  by  tho  magnitude  and  phase  of  the  reflected  indirect  ray. 
Reflectivity  appears  in  equation  9  as  the  product  of  roughness  and 
electromagnetic  reflection  parameters 

Tan  =  Cs  1  e_J° 


12. 


Fig.  5.  Pattern  propagation  factor  for  shipborne 
U.H.F.  radar. 
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Fig.  6.  Pattern  propagation  factor  for  shipborne  I  - 
band  radar. 
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SIGNAL 

POWER 

(dBm) 


RANGE  (n.m.) 


Fig.  7.  Signal  level  (dBm)  for  radar  of  fig.  5 
standard  atmospheric  conditions. 


14 


where  aa  =  HH  or  VV  represents  horizontal  or  vertical  polarization  and  0 
is  the  phase  difference  on  reflection.  In  the  program  a  semi-empirical 
expression  is  used  for  the  roughness  factor 


r  =  exp (-2  s1) 

=  exp(-1.2732  s) 


s<0. 6366  13a. 

s >0.6366  13b. 


where 


s  =  ji  siny/2X 


14. 


for  significant  wave  height  Equation  13  is  in  good  agreement  with 

experimental  data  over  the  range  of  grazing  angles  normally  of 
interest.  Figure  8  shows  the  difference  between  equations  13a/13b  used 
in  the  program  and  the  uniform  Gaussian  shape  (equation  13a)  used  in 
most  radar  models.  The  use  of  a  Gaussian  roughness  factor  leads  to  an 
underestimate  of  multipath  lobing  at  high  grazing  angles  and/or  sea 
states.  This  difference  can  be  significant  even  at  U1IF  frequencies. 


The  electromagnetic  part  of  the  reflection  coefficient  in 
equation  12  is  calculated  for  each  polarization  from  the  expressions 


(Pse"j%n  -- 

(pse  =■ 


siny  -  (ec  -  cos1y)0’, 
siny  t  Uc  -  cosJy)"'* 
t'csini  -  (cc  -  cosJy)0,5 
ecsir.y  +  (ec  -  cos*y)°‘5 


15a. 

15b. 


in  which  e„  is  the  complex  dielectric  constant  of  sea  water  at  the  air- 
c 

sea  interface: 


«c  -  «i  ‘  16, 

with  and  ej  being  the  real  and  imaginary  parts  respectively. 

Default  values  of  interfacial  sea  water  parameters  are  (i)  salinity  = 
3.4%  and  (ii)  temperature  »  sea  level  air  temperature.  Real  and 

imaginary  parts  of  ec  are  calculated  from  fits  of  experimental  data 

from  100  MHz  to  200  GHz.  Correct  allowance  for  temperature  and 

salinity  variations  is  more  important  for  vertical  polarization  where 
the  indirect  ray  grazing  angle  may  be  near  the  quasi-Brewster  angle  for 
high-altitude  or  short-range  targets. 


(GHz)  *  H  *  sin  (y  ) 


roughness  factor  (r)  for  pure  gaussian 
ssian /exponential  models. 
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4  NON-STANDARD  ENVIRONMENTAL  CONDITIONS 

In  the  absence  of  ducting  or  precipitation,  the  environmental 
parameters  which  have  the  greatest  effect  on  radar  performance  are 
humidity  and  sea  state. 

4.1  Humidity 

Absolute  water  vapour  content  is  calculated  in  the  program  from 
temporature  and  relative  humidity,  using  an  empi.  cal  form  of  the 
Clausius-Clapeyron  relation  (ref  2).  A  plot  of  attenuation  at  various 
humidities  is  given  in  figure  3  for  frequencies  from  100  MHz  to  25  GHz. 
The  effect  of  humidity  is  significant  only  above  1  GHz.  At  lower 
frequencies  the  molecular  absorption  is  dominated  by  oxygen  which  has 
negligible  variability  at  sen  level.  Absolute  humidity  and  atmospheric 
prossuro  are  also  used,  together  with  molecular  polarizations  of 
oxygen,  nitrogon  and  water,  to  calculate  refractivity  and  hence  the 
effective  oarth's  radius  as  described  in  reference  2. 

Figures  9  and  10  show  the  effect  of  humidity  variations  on  tho 
performance  of  a  hypothetical  5  GHz  surface  search  radar. At  high 
humidity  the  radar  horizon  is  increased  by  about  10%  as  is  the  poak  of 
the  first  lobe  for  a  1000  ft  target.  The  50%  increase  in  antenna 
temperature  does  not  show  up  in  the  noise  level  since  a  high  receiver 
noise  figure  (10  dB)  was  used.  An  additional  1  dB  increase  in  absorption 
near  the  horizon  is  also  not  apparent  due  to  the  simultaneous  shift  in 
the  lobe  structure. 

4.2  Sea  State 

Sea  state  contributes  (i)  as  a  clutter  return  and  (ii)  to  the 
multipath  lobing  through  the  roughness  factor  in  equation  12.  The 
program  assumes  a  simple  relation  between  significant  wave  height,  in 
feet,  and  sea  state 


H1/3  =  0.5  SS*  17. 

where  SS  is  the  Douglas  sea  state.  Figure  20  indicates  that  this  is 
reasonable  for  sea  states  up  to  7.  The  program  also  accepts  non¬ 
integer  sea  states  to  facilitate  the  investigation  of  smaller 
variations  in  wave  height. 
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5000  MHz 
250  kW 
10  dB 
1  MICROSEC 
1  MICROSEC 
50  FEET 


ANTENNA  TILT 


2  DEGREES 


HORIZONTAL  BEAM  WIDTH  r  1  DEGREE 

VERTICAL  BEAM  WIDTH  =  10  DEGREES 


VERTICAL  SIDELOBE  LEVEL  =  17-5  dB 


ANTENNA  NOISE  TEMPERATURE  =  108  KELVIN 

POWER  GAIN  =  35  dB 

TARGET  RADAR  CROSS  SECTION  =  1  METRE A  2 


DIELECTRIC  CONSTANT  OF  SEA 
(REAL  PART) 

(IMAGINARY  PART) 
RECEIVING  LINE  LOSSES 
POLARIZATION 

SEA  STATE  (DOUGLAS  SCALE) 
TEMPERATURE  AT  SEA  LEVEL 
SEA  LEVEL  HUMIDITY 
ATMOSPHERIC  ATTENUATION 


=  64-347 
:  36-146 

=  4  dB 

:  '  h 
=  3 

=  30  DEGREES  CELSIUS 

;  80% 

z  0-029893  dB/n.m. 


Fig.  9.  5GHz  radar  return  for  1m2  target  at  1000ft. 
T=30°,  humidity  =80%  and  state  3. 
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RECEIVER  NOISE  FIGURE 
COMPRESSED  PULSE  WIDTH 
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TARGET  RADAR  CROSS  SECTION 
DIELECTRIC  CONSTANT  OF  SEA 
(REAL  PART) 

(IMAGINARY  PART) 
RECEIVING  LINE  LOSSES 
POLARIZATION 

SEA  STATE  (DOUGLAS  SCALE) 
TEMPERATURE  AT  SEA  LEVE 
SEA  LEVEL  HUMIDITY 
ATMOSPHERIC  ATTENUATION 


=  5000  MHz 
=  2 SO  kW 
=  10  dB 
=  1  MICROSEC 
=  1  MICROSEC 
=  50  FEET 
=  2  DEGREES 
=  1  DEGREE 
=  10  DEGREES 
=  17-SdB 
=  71  KELVIN 
=  35  dB 
=  1  METRE  ~  2 

=  64*347 

=  36*146 
:  4  dB 
=  h 
=  3 

=  30  DEGREES  CELSIUS 
=  10°/. 

=  0*017670  dB/nm. 


Fig.  10.  5 GHz  radar  return  for  1m2  target  at  1000ft. 
T  =  30°  humidity  =10%  and  sea  state  3. 
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Figures  11  and  12  show  the  effect  of  sea  state  on  the  radar 
described  above.  Standard  conditions  are  assumed  with  sea  states  0  and 
5  respectively.  The  clutter  routine  used  is  based  on  a  semi-empirical 
fit  to  Nathanson's  data  (ref  6)  averaged  over  wind  direction,  as  shown 
in  the  sample  plot  of  figure  19.  A  small  return  at  sea  state  zero  is  a 
result  of  the  minimum  wave  height  (0.25  ft)  considered  in  the 
subroutine,  which  also  agrees  with  Nathanson's  data  for  sea  state  zero. 

For  target  ranges  between  5  and  30  n. miles  the  main  effect  of 
the  increased  sea  state  is  to  reduce  the  multipath  reflectivity.  This 
corresponds  to  the  region  of  maximum  slope  for  the  roughness  factor  in 
figure  8  and  results  in  an  8  dB  reduction  in  the  signal  level  for  the 
lobe  maxima.  The  practical  result  is  that  the  detection  range  for  a  1 
square  metre,  1000  ft  target  varies  between  25  and  50  n. miles  with 
normal  variations  in  environmental  parameters  and  system  degradation. 

At  UIIF  frequencies  the  sen-state  dependence  of  signal  return  is 
importart  only  for  higher  elevation  angles.  Figures  13  and  14  show  the 
variations  in  signal  return  for  a  sphere  droped  from  20,000  ft  at  10 
n. miles  ground  range,  which  is  a  typical  VCD  calibration  profile. 
Above  4,000  ft  the  low  sea  state  returns  are  5-10  dB  higher.  A 
Gaussian  model  for  the  roughness  factor  (r)  predicts  1-2  dB  less  lobing 
(figure  15)  for  high  sea  states  than  the  model  used  here.  Figure  8 
suggests  that  at  higher  frequencies  this  would  be  about  3  dB. 

5  PROBABILITY  OF  DETECTION 

A  second  program  uses  the  calculated  signal,  clutter  and  noise 
power  levels  to  compute  probability  of  detection  for  a  single  scan  with 
either  operator-  or  auto-detection.  This  is  the  same  as  the  blip-scan 
ratio  with  unity  operator  efficiency  if  the  integration  time  is  less 
than  the  scan  time  of  the  antenna,  but  greater  than  the  time  required 
to  scan  one  beamwidth.  The  main  routines  are  listed  in  Annex  II  and 
described  in  detail  in  reference  2. 

Target  fluctuation  is  synthesized  by  using  a  generalized  chi- 
squared  function  with  distribution  parameter  K,  which  may  take  any 
value  greater  than  zero.  Non-fluctuating  targets  correspond  to  K>>N 
where  N  is  the  number  of  pulses  integrated.  Weinstock  cases  are 
incorported  using  0.3<K<2.0,  and  Swerling  cases  I-IV  using  K=1,N,2,2N 
respectively. 
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Fig.  11.  5Gri2.  radar  return  for  1m2  target  at  1000  ft. 
T  =  15°,  humidity  =  60%  and  sea  state  0. 


FREQUENCY 
peak  POWER 
RECEIVER  NOISE  FIGURE 
COMPRESSED  pulse  width 
radiated  pulse  width 
antenna  height 

ANTENNA  TILT 
HORIZONTAL  BEAM  WIDTH 
'VERTICAL  BEAM  WIDTH 
ANTENNA  NOISE  TEMPERATURE 
TARGET  RADAR  CROSS  SECTION 

Dielectric  constant  of  sea 

(REAL  PART) 

(IMAGINARY  PART) 
RECEIVING  LINE  LOSSES 
POLARIZATION 

SEA  STATE  (DOUGLAS  SCALE) 
TEMPERATURE  AT  SEA  LEVEL 
SEA  LEVEL  HUMIDITY 
ATMOSPHERIC  ATTENUATION 


=  5000  MH2 
=  250  kW 
=  10  dB 
=  1  MICROSEC 
=  1  MICROSEC 
=  50  FEET 
=  2  DEGREES 
=  1  DEGREE 
=  10  DEGREES 
-  87  KELVIN 
=  1  METRE  ^2 

66*644 

35-883 

4  dB 
h 

5 

15  DEGREES  CELSIUS 
60% 

0*02267  dB/am. 


Fig.  12.  5  GHz  radar  return  for  1m2  target  at  1000 
T-  15  ,  humidity  =  60%  and  sea  state  5 
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Fig.  13.  UHF  radar  return  for  vertical  sphere  drop 
at  10  n.  miles.  Relative  drift  =  1  n.  mile  per 
40  K  feet,  sea  state  5. 
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Fig.  14.  UHF  radar  return  for  vertical  sphere  drop, 
sea  state  0. 
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Fig.  15.  UHF  radar  return  for  vertical  sphere 
drop,  sea  state  5,  gaussian  roughness 
factor. 
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Figures  16  and  17  are  plots  of  the  output  of  this  program  for 
processed  S/N  ratios  -15  to  30  dB  and  typical  values  of  N'=10  pulses  and 
PFA=10-<’.  Figure  17  shows  the  results  for  the  slowly  fluctuating 
Weinstock  targets.  At  90%  probability  of  detection  there  is  10-20  dB 
difference  in  the  required  signal-to-noi se  ratio  amongst  these  typical 
target  types.  That  is  to  say,  the  measured  VCD  contour  refers  only  to 
the  class  of  target  used  in  the  VCD  calibration.  The  difference  is 
considerably  less  i  f  R50  is  used. 

A  sample  of  the  graphical  output  of  program  II  is  given  in 
figure  18  for  the  5  GHz  radar  of  figures  9-12,  under  standard 
atmospheric  conditions,  and  incorporating  the  additional  losses 
described  in  section  2.4.  The  detection  probability  beyond  the  clutter 
horizon  (9  n.milos)  is  dominated  by  receiver  noise  considerations, 
while  there  is  nearly  100%  probability  of  paint  at  shorter  ranges  duo 
to  the  high  signal-to-cluttcr  ratio. 

Program  II  assumes  Rayleigh  statistics  for  surface  clutter. 
Over  the  range  in  which  this  is  a  valid  assumption  (r  >  0.25  pSec,  Ojj  > 

1  degree),  the  clutter  subroutines  give  +  5  dB  agreement  which  is  of 
the  order  of  the  source  data  accuracy.  This  level  of  agreement  is 
acceptable  considering  the  wide  range  of  clutter  RCS  encountered  in 
routine  calculations  (figure  19). 
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Fig.  16.  Blip-scan  ratio  for  fluctuating  targets. 
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chi-square  case  (K=3). 
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Fig.  17  Blip- scan  ratio  for 
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Fig.  18.  Blip -scan  ratio  using  radar  parameters  from 
figures  9-12,  under  standard  atmospheric 
conditions,  sea  state  3. 
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Wind  Waves  at  Sea 
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Figure  20:  Relationships  between  environmental  parameters  for  waves 
generated  by  local  wind.  (Source:  reference  6) 
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ANNEX  I:  Antenna  Pattern  Functions 

The  required  inputs  for  the  antenna  pattern  function  routines  are 
i)  3dB  beamwidths,  ii)  level  of  first  sidelobe  and  iii)  pencil  or 
cosecant-squared  pattern.  Pencil  patterns  are  calculated  according  to  the 
normalized  functions 

f(6)  =  sin  n(u2-B; ) 0  * 5  1-1 

f0n(u2-B*)“-* 

with 

f  _  sinh  tiB 
0  “  nB 

or  equivalently 

f(Q)  =  sinh  niB^-u2)0'*  1-2 

f0n(B*-u*) 0 • 5 

where  u=d.sin0/A  and  B  is  related  to  the  sidelobe  voltage  ratio  q 

q  =  4.603-si-n-h_nB_  1-3 

nB 

B»0  corresponds  to  a  sin  x/x  pattern  with  13.26  dB  sidelobe.  The  BASIC 
version  of  the  program  uses  a  cubic  fit  of  tabulated  solutions  for  the  above 
relation  and  for  the  aperture  efficiency.  The  numerical  accuracy  is 
comparable  with  that  of  the  recursive  algorithm  described  in  reference  2. 
Aperture  blockage  is  not  calculated  explicitly  but  can  be  approximated  using 
observed  sidelobe  levels  and  a  cosec  weighting  factor  as  plotted  in  figure 
1-1.  This  is  equivalent  to  assuming  that  the  envelope  decays  as  1/u  after  the 
obstacle  contribution  is  added.  Since  simple  parametric  forms  of  f(0)  are 
generally  adequate  only  for  the  first  few  sidelobes,  this  treatment  should  be 
reasonable  for  a  first-order  correction. 

The  program  is  easily  modified  to  use  a  lambda  function  subroutine 
for  f(0)  .  This  form  is  appropriate  for  an  aperture  with  distribution 
[l-(r/r0) 1]a“1  so  that  f(u)  is  a  lambda  function  of  the  first  kind 
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f(u)  =  2ap(l+a)-a—  -  1-4 

ua 

where  Ja(u)  is  the  Bessel  function  of  order  a.  The  subroutine  uses  a 
truncation  of  the  infinite  series  for  Ja(u).  Figure  1-2  shows  the 
output  of  this  subroutine,  with  the  series  truncated  at  u100  for  a=  0.5 
(sin  x/x) ,  a=  1.0  (uniformly  illuminated  circular  aperture)  and  a=  2.0. 
Circular  apertures  typically  have  a>1.0  (V>17.6  dB) .  Apart  from 

computational  simplicity,  tho  reason  that  the  modified  sin  u/u  pattern 
is  the  default  form  for  £(0)  is  illustrated  in  figure  1-3  where  the 
infinite  series  has  been  truncated  at  u^n  for  n  =  10,15,20  and  50. 
There  is  negligible  difference  in  the  position  of  tho  first  sidelobe 
for  the  two  expressions.  Tho  main  difference  is  that  the  sidolobe 
return  for  the  Bessel  function  form  diverges  rapidly  beyond  the  first 
sidolobe  unless  a  large  number  of  terms  is  included.  This  is  not 

convenient  for  routine  calculations. 


ANGLE  OFF  BORESIGHT - *- 

(degrees) 


Fig.  1-1.  Modified  sin  u/u  antenna  pattern  function 
for  9  =  1°,  B  =  0. 


0 


ANGLE  OFF  BORESIGHT  (degrees)  - ► 


Fig.  1-3.  Normalized  Bessel  function  (Ji(u)/u)  with 
series  truncated  at  u2n  with  n  =10,  15,  20 
and  50.  Modified  sin  u/u  with  17*6  dB 
sidelobe  shown  for  comparison.  BW  =  1°, 
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Annex  II: 

Probability  of  Detection  Subroutines 
for  Generalized  Chi-squared  Target 


1180  REM  I/P  to  1st  subroutine  (1320)  is  PFA  (P0), number  of  pulses  int'd  (N) 
1190  REM  1st  subroutine  is  used  to  calculate  threshold  for  fixed  threshold  det'r 
1195  REM  MMMMMMMMMtMMMMmMMMWMMItMMMMMMm 
1210  Y1=EXP(-Y0) 

1220  Y2=Y1 
1230  NO  =  1 

1 240  IF  N<=1  THEN  1290 

1250  Y1-Y1*Y0/N0 

1260  Y2=Y2+Y1 

1270  N0=N0+1 

1280  IF  N0<N  THEN  1250 

1290  RETURN 

1300  REM 

1310  REM 

1315  rem 

1320  REM  Start  of  first  routine  which  calculates  threshold  (YO)  from  PFA  (P0) 
1330  Y0=-1*LOG(P0) 

13*10  IF  N<=1  THEN  11110 

1350  REM  First  estimate  for  Y0 

1360  Y0=(SQR(-L0G( P0) )+SQR(N)-1 ) *SQR(-L0G( P0) )+N-SQR(N) 

1370  GOSUB  1210 

1380  D0=LOG( Y2/P0) *Y2/Y1 

1390  Y0=Y0+D0 

1H00  IF  ABS(D0/Y0)=>3.0E-7  THEN  1370 

1l)10  REM  Y0  is  threshold  for  fixed  threshold  detector 

11)20  RETURN 

11)25  REM  VMmmiimOfiMmtiiMOOMMIlllllllOllllMtOMMtMIHUKHHUHHHMiHHUHiO 
11)30  REM  Start  of  routine  to  calculate  Pd  from  signal- to-noise  (X) 

11)1)0  GOSUB  1210 
11)50  X=10A(X/10) 

11)60  SI  =Y2 

1D70  X1=(K/(K+N»X))AK 

11)80  X2=X1 

ID 90  X3=X*N/(K+N*X) 

1500  D1=K-1 

1510  M=1 

1520  R1 =0 

1530  L1=R1 

15*10  Y1=Y1/N0*Y0 

1550  Y2=Y2+Y1 

1560  S1=S1+Y1*( 1-X2) 

1570  El =( 1-X2) *( 1-Y2) 

1580  R1=S1+E1 

1590  X1=(D1+M)/M*X1*X3 

1600  X2=X2+X1 

1610  H=M+1 

1620  N0=N0+1 

1630  IF  ABS( 1-L1/R1 )=>3*OE-7  THEN  1530 
161)0  IF  El  =>0.1  THEN  1530 
1650  REM  rl  is  the  blip-scan  ratio 
1660  RETURN 

1700  rem 
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ANNEX  III:  Solutions  to  Multipath  Geometry 

The  pattern  propagation  factor  of  equation  8  is  evaluated  a 
points  in  the  interference  region  according  to  the  expression 

F  =  ffe^ltl  +  x2+  2x  cos  (0  +  2n6/X)  ] 0  ‘ 5  I 

where  0  is  the  phase  difference  on  reflection,  >.  is  the  wavelength 
8  is  the  path  difference  for  the  direct  and  indirect  rays: 

4R,R^  sin2y 
5  -  — L-A - 

Rf  +  R2  +  R 

where  R,  R^  and  R2  are  as  shown  in  figure  1  of  the  main  body  of 
note.  To  determine  these  slant  ranges,  the  corresponding  ground 
is  first  evaluated 

Gx  =  G/2  -  p  sin()j/3) 


with 


P  «  /  Whj)  +  (G/2) 2 


5  =  sin--1 


2aeG(h2-h1) 

1>5 


The  slant  range  is  then 


R  =  /oio-h^)  *  +4 ( a )  (ae+h2)  sin2 1  (Gj+Gj) /2ae ] 

Rj  (R2)  is  obtained  by  substituting  h2=0  (h2=0)  and  G2=0  (Gj=0) . 
solutions  for  ,  92  and  y  are: 

_i  2ae(h2-h1)+h2-hJ-R2 
G2  =  sin  1  e  1  1  11 

2(ae+h1)R 


G 


2  = 


sin 


-1 


2aeRl+Rl+Rl 

2(ae+h1)R1 


sin 


-1  2aehl+hl-Rl 
2aeRl 


t  all 

III-l 

and 

III-2 

this 

range 

III-3 

III-4 

I II- 5 

III-6 

The 

1 1 1-7 

III-8 


III-9 
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